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Abstract 

The single-particle inclusive differential cross-section for a reaction a + 
b — > c + X is written as the imaginary part of a correlation function in 
a forward scattering amplitude for a + b ^ a + b in a modified effective 
theory. In this modified theory the interaction Hamiltonian Hj equals 
Hj in the original theory up to a certain time. Then there is a sign 
change and Hi becomes nonlocal. This is worked out in detail for scalar 
field models and for QED plus the abelian gluon model. A suitable path 
integral for direct calculations of inclusive cross sections is presented. 
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1 Introduction 



In this article we consider inclusive cross sections, i.e. reactions of the type 

a( Pl ) + b(p 2 )^c(p 3 )+X, (1.1) 

where a, b, c are particles and X stands for the unobserved remaining reaction prod- 
ucts. We will present a general method which allows us to write the inclusive dif- 
ferential cross section p®da(a + b — > c + X)/d 3 p 3 as imaginary part of either a 
current-current or a field-field correlation function in a forward scattering ampli- 
tude a + b —>■ a + b in a modified theory. Let 

H = H + Hj (1.2) 

be the Hamiltonian of the original theory, with H and Hj the free and interaction 
parts, respectively. Then the modified theory is described by 

H = H + H I} (1.3) 

where Hj, obtained from Hj in a well defined procedure, is discontinuous in time 
and nonlocal in space. The modified theory is constructed in such a way that both 
its incoming and outgoing states are equal to the incoming states of the original 
theory. Thus the S-matrix of the modified theory equals the unit operator. Here and 
in the following we always work in the Heisenberg picture of quantum mechanics. 

Our article is organised as follows. In section 2 we recall the basic relations for 
single inclusive cross sections. In section 3 we present our general formalism for the 
modified effective theory in the case of scalar fields. Quantum electrodynamics with 
massive photons and the theory of quarks interacting with abelian gluons are con- 
sidered in section 4. We discuss some properties of the modified Hamilton operator 
and derive a path integral representation for inclusive cross sections in the abelian 
gluon model. These techniques are then applied to the cross section e + + e~ — » q + X 
as a specific example. We compare our techniques with the Schwinger-Keldysh for- 
malism |IJ, described e.g. in 0, and with Mueller's treatment [[3] of inclusive cross 
sections using the generalised optical theorem for 3^3 scattering in section 5 which 
contains also our conclusions. 



2 Single Inclusive Cross Sections 

In this section we recall some basic relations for inclusive cross sections. Our notation 
follows [|, [|. Let us consider a single-particle inclusive reaction, i.e 

a( Pl ) + b(p 2 )^c(p 3 )+X(p x ). (2.1) 
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To take a simple case, let a, b, c be spinless particles with masses m a} rrib } m c . The 
cm. energy squared is s = (pi + p 2 ) 2 - We assume c(p 3 ) to differ in type and/or 
momentum from a(pi) and b(p2). We use the covariant normalisation for our state 
vectors 

( a(p[) | o(p0 ) = (2vr) 3 • 2rf • 6®(jZ' -Pi) (2.2) 
and similarly for b, c. The S'-matrix element for the reaction ( |2.1| ) is given as 



S fi = ( c(p 3 ),X(p x ),out \ a(p 1 ),b(p 2 ),in) 

= iZ: 1 ' 2 j d A x e^ 3 ( X( Px ),out | j c (x) | a( Pl ),b(p 2 ),in ). (2.3) 

Here we have applied the reduction formula for particle c in the final state. Let 
4> c { x ) be a suitable interpolating field for c and Z c the corresponding wave function 
renormalisation constant. The current j c {x) is defined as 

Ux) = (a x + m 2 )(j) c (x). (2.4) 

The T-matrix element is obtained from the S'-matrix element via 

Sfi = 5 fi + i(2ir) 4 5 ( - A \p 1 +p2-p3-Px)Tfi, 

T fi = Z-V 2 { X( Px ),out | j c (0) | a(p 1 ),b(p 2 ),in ). (2.5) 

The single-particle inclusive cross section fmcips) is defined by 

d 3 a 

findPs) ■■= p° 3 -^-(a + b c + X) 

^ 4(2,)^( 1 s ,<, ma ¥ (2r)4f ' 4,(; " +K " ; ' 3 "^ )|r/ ' 12 ' (2 ' 6) 

w(x, y, z) = [x 2 + y 2 + z 2 - 2xy - 2xz - 2yz} 1/2 . (2.7) 

In the usual way the sum over all states \X, out ) in ( |2 .6|) can be carried out using 
completeness and translational invariance 
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4(2ir) 3 w(s, m\, ml) 
Z- l J d A x e ipsx (a( Pl ),b(P2),in \jl(0)j c (x)\a( Pl ), b{p 2 ), in > (2.8) 

= 2( 27 r)3 w (Ug,^) ImC ^^), (2.9) 

C(pi,p 2 ,p 3 ) = i J d A x e ipsx M(x), (2.10) 

M(x) = Z; 1 (a(p l ),b(p 2 ),tn\jl(0)j c (x)\a(p l ),b(p 2 ), m }e(-x°). (2.11) 

Here 9 (z) is the usual step function. 



2 



Alternatively we can write 

M(x) = Z^(n y + ml) (a z + m 2 c ) 

(a(px),b(p 2 ),in | (f)l(y)(j) c (x + z)\a( Pl ) } b(p 2 ),in )9(-x°)\ 

I y — >0 — , z^0~ 

(2.12) 

Here the limit y — > 0~, z — > 0~ is to be understood as follows: We first require y° < 
and z° < and perform the differentiations with respect to y and z. Afterwards we 
take the limit y — > and z — > 0. 

The amplitude C{pi,p2,ps) will play a central role in the following and we will 
be able to write it as a current- current, respectively a field-field correlation function 
in a forward scattering amplitude, but in a certain modified effective theory. 



3 Modified Effective Theory for Scalar Fields 

Let us assume that the basic dynamical variables of the original theory are the oper- 
ators for unrenormalised scalar fields <fii(x) and their conjugate canonical momenta 
Ilj(x) (i = 1,...,N). For simplicity we assume that EE* (or) = <f>i(x) holds. We denote 
4>i(x), IIj(x) collectively as Let H be the Hamiltonian of the system which we 

split into a free part H and an interaction part Hj which may depend explicitly on 
the time t, but should not involve time derivatives of IL(x) 

H(t, $(£ , t)) = #,,($(£, t)) + Hi(t, $(f , t)). (3.1) 

Besides the interacting fields and momenta $ free fields and momenta $^°- ) are 
considered with the corresponding Hamiltonian H . Here the mass parameters in 
H are taken to be the ones of the asymptotic particles. 

We assume now as usual (cf. e.g. 0, 0) that there existQ unitary operators U(t) 
that realize the time- dependent canonical transformations relating $ to $^°- ) 

$(f,t) = 17 _1 (*) ^ (0) t) I7(t) . (3.2) 

Taking as boundary condition 

$(f,0) = $ (0) (:r,0) (3.3) 



we get 



d t U(t) = -iHj(t,^°\x,t))U(t), 
U(0) = 1, (3.4) 



2 Of course this is only true in the regularised theory, i.e. for finite ultraviolet cutoff. 
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= £(-*)" / dt 



n=0 



dU 



dt n H I {t l )...H I {t r , 



We define furthermore 



U(t 2 ,t 1 ) = U(t 2 )U-\t 1 ). 



For t 2 > ti we have 

U(t 2 ,ti) = Texp 



— i 



t-2 



(3.5) 
(3.6) 

(3.7) 



where T means time-ordering. 

Let us recall the LSZ formalism [|, ||, [7|]. Assuming for simplicity the particles 
a and b to carry the quantum numbers of some fundamental hermitian scalar fields 
<p a , <pb, we define operators 

A( Pl ,x°) = iZ- 1 ' 2 j d*x e^ x Z°Ux) 

= iZ~ 1 / 2 J d 3 x e^(n a (x) - ip\4>a{x)\ 

B{p 2 ,x°) = iZ~ l/2 J d 3 x e ip2X K°<Pb{x) 

= iZl xl2 J d 3 x e ip2X {U b {x) - ip° 2 (j) b (x)), (3.8) 

where Z a b are the wave function renormalisation constants. The LSZ formalism has 
as basic assumption that for t — > ±oo the (hermitian conjugates of the) operators 
of ( |3.8j ) converge in the weak sense to the annihilation (creation) operators of out 
and in-states, respectively, for instance 



lim A\ Pl M)B\p 2 M) I ) = | a( Pl ),b(p 2 )^ ). 

£1,2— »±0O 



(3.9) 



Now we return to the single inclusive cross section (|2.8| ), where we have to cal- 
culate the matrix element ( |2.11| ) 

M(x) = Z~ l { a(pi),b(p 2 ),in \ jl(0)j c (x) | a(pi) , b{p 2 ) , in )9(-x°) 
= Z- 1 lim ( | B(p 2 , -t' 2 )A{p x , -t' x ) 

jl(0)j c (x)Ai(pi, -ti)B\p 2 , -t 2 ) | )9{-x°). (3.10) 

We assume the current j c {x) to be expressible in terms of the fields and their con- 
jugate momenta but not involving their time derivatives. Then we have 

U(t)j c (t,^x,t))U-\t) =j c (t,&°\x,t)) =:ji°\x). (3.11) 
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Furthermore we have from ( |3.2| ) and 

U(t)A( Pl ,t-^(x,t))U-\t) = A( Pl ,t;&°\x,t)) =: A<°>(t) (3.12) 

and the same for B(p2,t). We note that the <3? ( - * ) satisfy the free field equations. In 
the case of a free scalar field with mass m they read 

<p {0 \x,t) = n(°)(x,t), 

fl (0) (x,t) = (A x -m 2 )^°\x,t). (3.13) 

This implies 

i (0) (£)=0, 5 (0) (t)=0. (3.14) 
Now Ai(x) (|3.10|) may be written as 

M{x) = lim Z~ x { | U-\-V) 

U{—T', -t' 2 )B^(-t' 2 )U(-t' 2 , -t'JA^i-t',) 

A t(0) (-ti)f/(-ti, -t 2 )i? t(0) H2M-i 2 , — T) 

U(—T) | )0(-a; Q ). (3.15) 

Here we have introduced further times T, T' and we assume without loss of generality 

T>t 2 > h, 

T'>t' 2 >t[. (3.16) 

With the usual assumption that the interaction is switched off adiabatically for 
t — > ±oo we get from the adiabatic theorem 

lim U(-T) | ) = e iLp - \ ), (3.17) 

T — >oo 

where |0) is the vacuum state of the free theory. Inserting everything in ( |3.15| ), we 

get 

M.(x) = lim lim Z^ 1 

ti,t'.-KxT,T'-Hx 6 

( o | U ^ T \-t' 2 )B^\-t' 2 )U{-t^-t\)A^\-t' l ) 
C/(-t / 1 ,0)itW(0)C/(0,x°)if (xMx ,-^ 
A^°\—ti)U(—ti, —t2)B^ \—t 2 )U(—t2, — T)| )9(-x°) 
[( | U{-T')U-\-T) | J]" 1 . (3.18) 

Following the time-dependence in Ai(x) from the right to the left, we start at time 
—T — > — oo and pass through operators of increasing time arguments until time 0. 
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Then the time sequence changes and we go back in time to time — T' — > — oo. In a 
usual matrix element the time arguments should increase instead. 

We will now show that we can write the matrix element M.(x) ( 3.18|) in the 



usual form, with time increasing from right to left, if we pass to operators U of the 
form ( |3.5| ) but with a modified interaction Hamiltonian Hj. As the time-sequence 
in M.(x) is correct up to t = 0, we request 

U(t) = 6(-t)U(+t) + 0(t)U(-t). (3.19) 

This gives with ( |3.4|) 

d t U(t) = -iHj(t)U(t), (3.20) 

Hj(t) = ei-^H^^^x^-e^H^-t^^^-t)). (3.21) 

For t > our modified interaction Hamiltonian Hj(t) depends on the free fields and 
momenta at time (—t). But we know how to express the free fields and momenta 
at time (—t) by their values at time t using the free field equations of motion. For 
, n(°> satisfying ( |3TT3| ) we get 

0(°)(f )X °) = f d 3 y A(x-y;m 2 )^o0 (o) (y) 

J y°=const. 



( d 3 y {A(x - y;m 2 )U {0 \y) + A(x - y; 

Jy°=const. '- 



m 2 



I y°=const. 

n(°)(f,x°) = (o) (f,x°) (3.22) 
with the usual commutator function for scalar fields of mass m 

A(z;m 2 ) = % f d R k(e- ikz -e ikz ), 



d R k = -^-8(k°)5(k 2 -m 2 ). (3.23) 

[ATI) 6 

Setting y° = t and x° = —t in (|3.22Q we get expressions for (f)(°> (x, —t) and n^ ) (x, —t) 
as linear functions of (j)^ '(y,t) and US°'{y, t). Thus, we can consider Hi(t) (|3.21| ) as 
a nonlocal functional of the dynamical variables t) at the same time t also 

for t > 

H I (t) = H I (t,&°\x,t)). (3.24) 
Using in addition (|3.14|) we can rewrite Ai(x) ( |3.18|) as follows 
M.(x) = lim lim Z7 1 

ti,t'.-+ooT,T'->oo 

(0\U(T'X)B {0 \t' 2 )U(t' 2 ,t[)A(°\t^ 

U^O^m^x^^Uix^-h) 
At(°)(-t 1 )t/(-ti,-t2)5 t(0) (-i2)^(-t2,-T) I )9(-x°) 

[( | U(T',-T) | )]~\ (3.25) 
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U(t',t) = U{t')U~ \t) = Texp 



(t' > t). (3.26) 



In the limit T, T' — > oo we get from (|3.25| ) 



A* fx) 



lim Z c 1 ( | T exp 

tit'— >oo 



dt'Hiit') 



B^\t' 2 )A^\t[)j^\0)ji°\x)A^(~t 1 )B^(-t 2 ) | 



| Texp 



— i 



dt'H^t') 







-i 



(3.27) 



Clearly, we can consider (|3.27|) as matrix element of the standard type but in the 
modified theory governed by the total Hamiltonian 



H(t, $(£, t)) = H {${x, *)) + Hj(t, Hx, t)), 



(3.28) 



(3.29) 



M{x) = Z-\{ a( Pl ),b(p 2 ),out | jt(0)j c {x) | a(p 1 ),b( P2 ),in ))6(-x°). (3.30) 



Here we denote by (( )) matrix elements in the modified theory. Using 
also write M. as 



Mix) 



m. 



(( a{pi),b{p 2 ),out | <f>l(y)(f) c (x + z) \ a{p x ), b(p 2 ), in ))9(-x°) 



we can 



\y^0-, z->0- 

(3.31) 



The S-matrix of the modified theory equals unity, due to the fact that the in- 
coming and outgoing states are now identical 



S=l. 



(3.32) 



In ( p.28 ) - ( |3.31| ) we have a main result of our paper. The matrix element M. is 
written in the standard form with an in-state to the right and an out-state to the 
left. The prize we have to pay is that we have to use the modified theory where the 
Hamiltonian H(t) has a sudden variation at t = and is nonlocal for t > 0. On the 
other hand we can use ( [3.28D - (|3.31|) as starting point to write down a path integral 
representation for Ai in the standard way. Below in section 4 we will do this for 
QED and the abelian gluon model as an example. 
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4 Inclusive Production in QED and an Abelian 
Gluon Model 



4.1 QED and the Abelian Gluon Model 

In this section inclusive reactions are considered in QED coupled to an abelian gluon 
model. An example for such a reaction is 

e + {pi) +e-(p 2 ) -> q(pa)+X. (4.1) 

Starting point is the Lagrangian describing the interaction of electrons of mass m 
and charge — e with a massive photon of mass A - to avoid any infrared divergences - 
and of two quark flavours of equal mass M and electric charge eQ q with the photon 
and a massive abelian gluon of mass \i. As Lagrangian we choose 

£ = - -G^G^ - — (9„G") 2 + -filG.G^ 
4 2r) 2 M 

4 M 2£ 2 ^ 

+ V Q £ -m + e 4) ^ + 9 Q ? -Mo - e Q 9 /-2or 3 <?) g, (4.2) 

where G^ denotes the abelian gluon field and G^ v = d^G u — d v G^ its field strength 
tensor, the photon field and = d^A u — dyA^ its field strength tensor, ip the 
electron field and q the quark field. For the vector boson part of £ we have chosen 
here a form advocated by Stueckelberg (cf. e.g.[0), where terms —l/2r] (d^G fJ- ) 2 and 
— l/2£ (<9 M v4^) 2 have been added for the abelian gluon and the photon field. We have 
a quark field with two flavours and r 3 is the usual Pauli matrix. All quantities in 
( |1.2j ) are the unrenormalised ones. The corresponding renormalised quantities are 
denoted either without any subscript in the case of constants or with the subscript 
R in the case of fields. 



The Lagrangian ( f4.2| ) has two discrete symmetries of the charge conjugation 
type, Ca and Cq, which guarantee that no further coupling terms need to be added 
in ( |4.2| ) and that photons and the abelian gluons do not mix. All this is explained 
in appendix A, where we also discuss the renormalisation procedure in our model. 

The canonical momenta of the vector fields of the theory are listed in (|B.2|) of 



appendix B. For finite £ an d r/ the relation ( |B.2| ) can be solved for OqA^ and OqG^ 
and the canonical formalism can be applied. 



The Hamiltonian corresponding to (|4.2j ) can be obtained as usual by a Legendre 
transformation 

H = Ho^+Htiy ), 
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H (y°) = Jd 3 yH (y), 

Hriy ) = fd'yHriy), (4.3) 



where we choose the free T~Lo{y) and the interaction part T~ti(y) of the Hamiltonian 
density as follows 

n {y) = -|n GO n£ - ^n G TP G - s^n^ 

-^u A0 u° A - in Ai n^ - dA'K + n%A° + - 1 -\ 2 a^ 

(frdi -m) V - q Qf£ -M) g. (4.4) 

= + 1 -5i?g,g» + ^77n G0 nS + ^a 2 ^ + ^n A0 n° 

—Smipif) — 5Mqq 

-e ^M + q(e Q q fl+g T 3 a)q } (4.5) 

with 

; 2 2 2 

dfx = fi -/j, , 
S v = 7] - r] , 
SX 2 = A 2 -A 2 , 

5m = m — mo, 

5M = M-M . (4.6) 

Only the interaction term has to be altered in order to construct the Hamilton 
operator for the modified theory according to the rules specified in section 3 

H(y°) = Hoiy^+Hjiy ), 

Hj(y°) = f<Py[0{- 1 / i )Hi(y)-e( 1 / ) )H I (v)] 

with y=(-y°,y). (4.7) 

In Hi(y) the fields at time —y° have to be substituted by the fields at time y° via the 
time shift relations analogous to ( |3.22| ). For the vector fields and the fermion fields 



these are given by (|B.6|) , (|B.7| ) and (|B.11|) in appendix B. Two points should be 



stressed concerning this modified Hamilton operator: First, due to the jump in the 
interaction, the effective theory has no time translation invariance, so that no longer 
energy conservation holds at every vertex. Second, the effective interaction contains 
derivatives. Thus, if we use perturbation theory for the evaluation of amplitudes 
in the effective theory, we have to be careful about contact terms in time ordered 
products, especially for the vector fields (for a general discussion of such contact 
terms see |J). 
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4.2 Inclusive Quark Production in e + e -Annihilation 



We want to calculate the single-quark- inclusive cross section of reaction (|4.1|) . Per- 
forming all the steps as described for scalar fields in section 3 we can write this 
inclusive cross section as imaginary part of an amplitude C ( cf. ( |2.8| , |2.10] , 2.12|) ). 
For unpolarised e~ and e + in the initial state and summation over spin and flavour 
of the final state quark this reads 



fmM = W^ ImC(?1,?2 ' ?3) ' (4 ' 8) 



C(pi,P2,Ps) = t J d 4 x e ip3X M(x), (4.9) 
M(x) = -Z- l Y.u{p,W z -M) 



e~ > e 



(e + (pi),e (p 2 ),in\ Tq A (x + z)q A (y) \ e + (p 1 ),e (p 2 ),in)Q(-x ) 

- M)u(p z )\ 

(4.10) 

The transition to the effective theory gives 
Mix) = -Z^YHPsWz-M) 

£'« e + ( Pl ),e-(p 2 ) | Tq A (x + z)q A (y) \ e + (p 1 ),e~(p 2 )}}6(-x ) 

S 4- ,8 — 

- M)u{pz\ 

I j/ — >0 , z— >0~ 

(4.11) 

Here Yl means the average over the spin states, A is the quark flavour index over 
which we sum. In ( (4.10j ) and ( |4.11| ) we have for convenience inserted a T-product 



symbol which puts q and q into the order required for x° < in the limit y — > . 

In order to calculate M(x) we expand it in powers of the electromagnetic cou- 
pling constant e. The leading term in the connected part of M is of order e 4 , cor- 
responding to the diagrams of figure 1 which can be classified as e + e~-annihilation 
(a) and e + e~-scattering (b) ones. 

Using the LSZ reduction formalism we get for the connected part of the annihi- 
lation matrix element M a 

H'<(e + (pi),e~(p 2 )| Tg A (a; + z)q A {y) | e + (pi), e~(p 2 ))) = 

e> ' e 

= -e 2 /^((0| T Jd 4 x 3 {9(-x° 3 )A^x 3 )e i{pi+P2)xa - 6(x 3 )A^x~ 3 )e i{pi+P2)x ~ 3 } 
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Figure 1: Annihilation diagrams (a) and scattering diagrams (b) for the amplitude 
M.(x) in order e 4 . In the middle blob arbitrary gluonic interactions are allowed. All 
lines and vertices are in the effective theory. Time runs from right to left. 



q A (x + z)q A (y) 

Jd^iei-xDA^x^e-^ 1 ^^ 4 - e{xl)A u (£ A )e' i ^ 1+p ^} |0)), 

(4.12) 

where we have written the matrix element with the help of the lepton tensor 

s + ,s _ 

= P&2+P1P2 -g^ipm + m 2 ). (4.13) 
Expanding further we get to lowest order in Q q e 

^'((e + (pi),e'(p2)\ Tq A (x + z)q A (y) | e + (p 1 ), e~(p 2 ))) 

S e+> S e~ 



((0| T jd%[0{-xl)A»{x z y^+P*)** -e(x° 3 )A^x~ 3 )e^ +p ^ 

'd 4 x 1 [6{-x° 1 )q{x 1 )A'(x 1 )q{x 1 ) - fl^?)^)^)^^)] 

q A (x + z)q A (y) 
d 4 x 2 [e(-x° 2 )q(x 2 )AXx 2 )q(x 2 ) - e(x° 2 )q(x 2 )A[x 2 )q(x 2 )} 

Jd 4 x A [e(-x° 4 )A u (x 4 )e~^ +p ^ -e(x° A )A u (£ 4 )e-^ +p ^] |0» = 
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-e 4 Q^- 



1 



1 



A 2 — it s — A 2 + ie 
T yc/ 4 xi[0(-x?)e l(pi+P2)xi g(xi)7 M g(xi) - ^(x°)e i(pi+P2)£l g(xi) 7/ ,g(x 1 )] 

q A (x + z)q A {y) 

U A x 2 [9{-x 2 )e- i ^ + ^q{x 2 ) lv q{x 2 ) - e{x 2 )e- i(pi+P2)i2 q{x 2 ) lu q{^)]\0)) , 

(4.14) 

so that our matrix element Ai a (x) can be written in the following way 



M a (x) = Z-'e^QT 



X 2 + ie 



J2HPs)(^z-M)e(-x ){(0\Tq A (x + z) 
ld%[6(-x° 2 )e-^ +p ^q(x 2 ) lu q(x 2 ) 



9(^2)7^5(^2)] 



rf 4 x 1 ^(-x°)e i(pi+P2)ll g(xi)7 M g(x 1 ) - ^(x?)e i(pi+P2)5 ' 1 g(xi) 7/ ,g(xi)] 

g A (2/)|0))HS-MH^)i 

l»-»0-, 2-^0 

(4.15) 

in terms of quark 6-point-functions in the quark-gluon-theory. In the next subsection 
this will be represented by a path integral in the modified abelian gluon model. 
The electromagnetic interaction of the incoming fermions could be separated, the 
calculation is sketched in Appendix C. 

The same procedure can be applied to the scattering matrix element Ai b (x) 
corresponding to the diagram of figure lb. Of course this diagram cannot give a 
contribution to the inclusive cross section. In higher orders the e + e~-scattering type 
diagrams give e.g. the contributions of 2-photon annihilation processes (cf. figure 2) 
to the inclusive cross section ( |4. 1|) 

e + + e~ -> e + + e~ + q + X. (4.16) 

As a simple check of our theoretical manipulations let us finally calculate the 
matrix element A4 a (x) in lowest order in the quark-gluon-coupling g, i.e. for g = 0. 
Decomposing the quark 6-point-functions with Wick's theorem into quark 2-point- 
functions and substituting for these the perturbative propagators ( |B.14|) we get for 
Ai a two contributions according to the two different possible contractions 

1 2 



M a (x) 



2e*Q 2 F u 



\ 2 + ie 



6(-x°) 



7^-Mh^ 



+ e 



i{p+P\+P2)x 



lv M-M) lv ] u(p 3 ). (4.17) 
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Figure 2: The diagrams of 2-photon annihilation processes for inclusive quark pro- 
duction (|4.1|) are shown in (a), the corresponding diagrams for M. ( [4.1 1|) in (b). 
Time runs from right to left. 



Inserting this in Q4.8 ) we get 



2(2ii) 3 w(s,m 2 ,m 2 ) q s — X 2 + ie 




P=P1+P2~P3 



(4.18) 



This is, of course, the standard result, which one obtains in considering to lowest 
order the reaction 



for two quark flavours of charge Q q . 

4.3 Path Integral Representation 

In the last subsection the matrix element Ai(x) and therefore the one-quark-inclusive 
cross-section could be expressed in terms of quark 6-point-functions after separating 
the electromagnetic interaction by a perturbative calculation. As the coupling g 
is not assumed to be small we will now derive a representation of these quark 6- 
point-functions suitable for non-perturbative calculations. For this we consider the 
Hamiltonian path integral in our effective theory obtained from the abelian gluon 
model. 



e7 + e 



q + q 



(4.19) 
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Any Green's function of the theory can be written as 



Tq{x 1 )...q{x 2 )...G x {x 3 )...U G {x 4 )... |0» 

= Z- X J V(G,U G ,q,q) q(x 1 )...q{x 2 )...G\x 3 )...U G (x 4 )... 

x exp{* f dSj (!%)<%) - H{y))} (4.20) 



with 



J V(G, n G , q, q) exp{* j dSj (1%)%) - H{y))} (4.21) 



and 



n(y)$(y) = U G (y)G p (y) + n q (y)q(y) + ^y)n q (y) 

= n G (y)G p {y) + ^q{y) lo d q(y). (4.22) 

Here the classical fields G and Tl G and the Grassmann fields q and g have to be 
inserted into the part of H(y) (|4.7|) which describes the abelian gluon model. Fur- 
thermore for y° > the classical and Grassmann fields have to be time-shifted 



according to (|B.6|) and ( |B.11|) so that we get 



IL(y)$(y)-H(y) = 

l Vo U GO (y)U° G (y) + d^(y)n° G (y) + lu Gi (y)n G (y) - G iQ {y)tf G {y) 



-- GlJ (y)G^(y) + -X 2 G,(y)G^y) + q(y) ^ y -M j q(y) 

- 0(-y°) 9o q(y) r 3 G(y) q(y) 

+ %°) 9o q(y) r 3 G(y) q(y). (4.23) 

For the modified abelian gluon model the canonical momenta U G can be inte- 
grated out explicitely to obtain a Lagrangian path integral 

J V(G,U G ,q,q) exp{* JdSj (U(y)^(y) - H(y))} 

= J V(G, q, q) exp{* J dSj C(y)} (4.24) 

with an effective Lagrangian C(y) 
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+ q{y){\?y-M,)q{y) 

-0(-y )9oq(y)T 3 a(y)q(y) + e{y°)g a q{y)^Q{y)q{y) 
-0{y°)\gl j d% 5«(y° - z\) J d% S^(y° - zl)q{z~ x )r, 

{-(-A(ii - y; Ag) 7o + 3)[A(z a - y; r) \ 2 ) - A(ii - y; Ag)])g(£i) 

g(i 2 )(-A(i 1 - y; Ag) 7 ° + ^ ^°[A(ii - y; ^) - A(f x - y; Ajj)]) 

+(-A(ii - y; Ag) 7i + ^^[A(fx - y;^) - A(ii - y; Ag)])g(ii) 

q(z~ 2 )(-A(z\ - y; Ag) 7 4 + -U ^[A(ii - y; i^Ag) - A(ii - y; Ajj)])} 

r 3 g(f 2 ). (4.25) 

The straightforward calculation leads to a four fermion coupling term, which results 
from the time-shift relations of the bosonic fields. 



As the Hamiltonian path integral ( [4.21 ) is quadratic in q, q, the fermionic fields 



can be integrated out. With the Green's function S F of the quark in the modified 
gluon background 

(iP-M Q )S F {z 1} z 2 ; G)-gJ d A z G{z x - z)S F (z, z 2 ; G) = -6^( Zl - z 2 ), (4.26) 

the quark propagator is given as 

(( | Tq( Zl )q(z 2 ) | 0)) = (( \s F ( Zl , z 2 ; G))), (4.27) 

where the brackets (( )) on the right hand side denote the average over all gluon 
fields with the measure dictated by the path-integral 



({F(G,U G ))) = Z- 1 Jv{G,U G )F(G,U G ) 



x exp [ i J dSj (^ n G0 n^ + d^U° G + ^U Gl U G - G i0 U G 

~^GijG tj + -AoG^G 1 ')] 
x (detH (iP-M -g J d A zG{ .; z))])' (4.28) 
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with 



fll/iJlfe) = e(-y 1 )5^(y 1 ~y 2 )T 3 (^(y 1 ) 



j S( yi - ~z; M)t,Q(~z)S(~z - y 2 ; M) 7 °. 

(4.29) 

Inserting the Green's function ( |4.27| ) into the matrix element ( |4.15| ) we get the 
matrix elements and therefore the inclusive cross section expressed explicitely in 
terms of a Hamiltonian or Lagrangian path integral, respectively. These expressions 
should be a convenient starting point for applying non-perturbative methods to an 
evaluation of inclusive cross sections. It should be possible, for instance, to generalise 



the methods of [11 to the case of this effective theory here. 



5 Comparison to Other Techniques and Conclu- 
sions 

In this section we give a brief comparison of our technique using time-shifted fields 
with the formalisms of Schwinger and Keldysh |l| ||| and Mueller ||. Then we will 
draw our conclusions. 



To discuss the Schwinger- Keldysh formalism we go back to (3.15p where the time 



sequence runs forward from — T to and then back to — T". Following Schwinger 
and Keldysh one inserts [/-matrices from to a large positive time T" and back to 




U(0,T")U(T",0) = 1, (5.1) 



Ai(x) = lim lim Z r , 1 

( o | U{-T\-t' 2 )B^\-t' 2 )U{-t , 2 ,-t' l )A^{-t' l ) 

f/(-t;,o)f/(o,n^(r^o)jf)(o)f/(o,x )j( )(x)f/(x ,-t 1 ) 

A t(o) ( _ tl)[/( _ tl; _^ )s t(o) ( _^ )c/( _ £2) _ T) i o )e( _ x 0) 

[( I U(T',-T) | (5.2) 

In the complex time plane we go along a path starting at — T, going to T" and then 
back to — X" (figure 3). This method is best used introducing the two field formalism 
P], H . It plays an important role in thermal field theory in the real time formulation 

11 
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t-plane 



—T 


rpll 














-V 





Figure 3: Time contour in the Keldysh formalism. 

In Mueller's formalism one starts with the 3-3 amplitude for a + b + c — > a + b + c 

( 6(2/3), c(pg), out I a(pi),b(p2),c(p 3 ),in ) = 

= 6 fi + i(2n) 4 6 {4) (pi + p 2 + p 3 - p[ -Pi- p' 3 )T fi . (5.3) 

Then one identifies a particular contribution to Im7}j, i.e. a particular discontinuity 
of T fi (figure 4). 

c c 
a a 

Figure 4: In (a) the discontinuity of Tfi is shown which is related to the inclusive 
cross section a + 5 — > c + X in Mueller's formalism (b). 

We should emphasise that our approach is quite distinct from these. In our 
effective theory the time contour runs from —00 to +00. The scattering matrix 
S = 1 (cf. p.32| ). If we insert the matrix element M. calculated in the effective 
theory according to ( |3.31| ) into ( |2.10| ) we get an expression which looks almost like a 
3— >3 scattering amplitude but not quite so because of the function 9(— x°) in ( |3.31| ). 

To conclude: We have presented in this article a novel technique for treating inclu- 
sive reactions: the method of the time-shifted fields. We have explained this method 
for theories with scalar fields and with electrons and quarks interacting with photons 
and abelian gluons. The generalisation to QCD should be straightforward and will 
be dealt with in future work. We have written the inclusive cross section as imagi- 
nary part of an amplitude C for which we have given a path integral representation. 
We have applied our formalism to the reaction e + + e~ — ► q + X and checked that in 
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lowest nontrivial order we get the correct result. Our hope is that our path integral 
representation will lead to a genuinely non-perturbative evaluation of inclusive cross 
sections at high energies along the lines of |TT] . The method of the time shifted fields 
could also be useful for general studies of inclusive reactions in QFT. One can, for 
instance, think of generalising Wilson's operator product expansion (OPE) |H| to 
our effective theory. Then our methods should allow a straightforward description 
of inclusive production of hadrons h in e + e~ annihilation at high energies 

e+ + e~ -> h + X (5.4) 

in terms of a genuine OPE. A similar approach should be possible for fracture 
functions for hadron-hadron or virtual photon (7*)-proton scattering 

h x + h 2 -> h 3 + X 

7*+p _> p + X. (5.5) 

The last reaction is, of course, related to diffractive deep inelastic scattering as 
observed at HERA |T7|]. Finally, our methods can also be applied to the treatment 
of non-equilibrium processes in thermal field theory. 
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Appendix 

A Relations for QED and the Abelian Gluon 
Model 

A.l Symmetries 

The Lagrangian (|4.2| ) has the following two symmetries of the charge conjugation 
type 

eCq T (x), 
-A^x), 

G»(x); (A.l) 
ip(x),, 



C A : if>(x) 
q(x) 

G"{x) 
C G : Mx) 
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q(x) -> -eq(x), 
A"(x) -> 

-> -G"(z). (A.2) 
Here C is the usual charge conjugation matrix for Dirac fields (cf. e.g. |12| ) and 

(A.3) 

The standard charge conjugation transformation is 

C = C A C G . (A.4) 

The symmetries Ca and Cq forbid coupling terms of the types AG, A 3 , AG 2 , A 2 G, 
G 3 , AG 3 , A 3 G in our Lagrangian Q4.2| ) and require the quark mass term to be flavour 
diagonal. Coupling terms of type A A , A 2 G 2 , G 4 have mass dimension 4 and could be 
included in C But we can leave them out since the renormalisation procedure does 
not require them. The reason is as in QED: The vector bosons couple to conserved 
currents and thus the superficial degree of divergence of vertex functions of four 
vector bosons is less than (cf. e.g. 0). 

A.2 Renormalisation 

We outline here the main steps in the renormalisation program of our theory. 

Let us define the unrenormalised vacuum polarisation tensors LT^ (k) of the pho- 
ton and n.ffi(k) of the gluon 

= ie 2 U^(k), (A.5) 
^astrm^v = iglufXk). (A.6) 



Here the blobs mean the sum of the one-particle irreducible diagrams. In the stan- 
dard way one shows that 11^ and H^y are purely transverse and independent of 
the "gauge parameters" r} and £ - Thus we write 

n$ G) (k) = i-g,uk 2 + k,K)rt A > G \k 2 ). (A.7) 

' (A) 

The unrenormalised photon propagator A^ reads 
iA™(*) 



"" "' k 2 [l + e 2 TK A )(k2)]-\ 2 + ie 

W» + (Co[l + e 2 rt A \k 2 )} - 1) k ;%,. ]. (A.8) 
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The renormalised photon mass A is obtained from 

A 2 [l + e 2 IT (A) (A 2 )]-A 2 = 0. (A.9) 
The wave function renormalisation constant Z$ of the photon is given by 

zi A) = [l + ejjll^A 2 )]- 1 . (A.10) 

It follows that the renormalised mass A, the renormalised charge e and the renor- 
malised gauge parameter £ are given by 

A 2 = zi A) \l (A.ll) 
e 2 = Z^el, (A.12) 
£ = (zi A Y%- (A.13) 

With that we get for the renormalised field and propagator 

*k(x) = (4 A) )-^(x) t (A.14) 



M) (k) - Zl_ 



k 2 [l+e 2 U ( c\k 2 )]-X 2 + ie 

\9r + m + e 2 Ui A \k 2 )} - l) k2 _ k f x ; + J (A.15) 

n^ 2 ) = n^(^ 2 )-n (A) (A 2 ). (A.i6) 

Note that with the conventional choice of Z% in ( |A.10|) the residue of the pole at 
k 2 = A 2 in the transverse term of the renormalised propagator is finite but not 
normalised to 1. 

For the gluon relations completely analogous to ( |A.8 ) - QA.16 ) hold with the 
replacements A — > G, eo — > go, \o — > fio and £o —* Vo- 

The renormalisation of the electron and quark fields is standard. Let —iYs$> (p) be 
the unrenormalised electron self energy, then the unrenormalised electron propagator 
S'M is 

S'^ip)- 1 = tf-m - EW(p) + ie. (A.17) 

The renormalised electron mass m and the wave function renormalisation constant 
Z 2 are given by 

= m-mo-S(m), (A. 18) 

Z? = [l-^W(/)^J-l (A.19) 



-ra 
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For the renormalised electron field ipn and propagator Sp R we have 

i/> R (x) = (Z^)-h(x), (A.20) 

S^Up)' 1 = /-ro-E^ (?)+«, (A.21) 

E^(p) = #[EW(p)-SW(m)]-(#-l)(/-m). (A.22) 



Relations analogous to ( A. 17 ) - ( |A.22| ) hold for quarks in our theory with the 



replacements ift — > q, m — > M . Note that also the complete quark propagator is 
flavour diagonal due to the Ca and Gq symmetries ([A.1|) ,( [^2"D . 



This concludes our brief discussion of the renormalisation of the model. 



B Free Fields and Time-shift Relations 

B.l Massive Vector Field 

A free vector field of mass \x can be described by a Lagrangian density 

C = -\g^ u - ^(d^f + \v 2 G»G» (B.l) 
with the field strength tensor G^ v = dpG v — dyG^. With the canonical momenta IP 



the Hamiltonian density is given as 

n = -^n n° - injr - d t G l u° + ir^c + ^g^ - \^ 2 g^. (b.3) 

Using the commutator function A 

A(x; fi 2 ) = i J d R k{e- ikx - e ikx ) (B.4) 

with 

^*=(03W (1) (* a -A* 2 ), (B.5) 
G^ at the space-time x can be obtained from G M at time y° by means of 
Gp(x) = f d 3 y[A(x-y;i2 2 )G fl (y)-A(x-y;i2 2 )Tl fl (y) + 

J yO=const. 



+\d»d p [A{x - y- W 2 ) - A(x - y; y?)]W{y) 

+8,A{x - y; W 2 )G°(y) - g ,d p A(x - y; /i 2 )G%) ], 

(B.6) 
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where the time- derivative of has been expressed in terms of and II M with the 
help of the free equations of motion. The derivatives are taken with respect to the 
first argument. In the case of rj = 1 (|B.6|) reduces to 



Gp{x) = I d 3 y A{x-y-f%G^y). (B.7) 

J y°=const. 

Furthermore we need for the perturbative calculation the propagator 

(0| TG^GM 10) = / ^j^"^- (B.8) 

B.2 Dirac Field 

Starting from the Lagrangian density for a free Dirac field ip of mass m 

C = -m)^, (B.9) 
the corresponding Hamiltonian density is given as 

H= -$(W§i -rn)il>. (B.10) 
The mapping of ip from y to another space-time point x is accomplished by 

ip(x) — —i d 3 yS(x — y;m)'y il)(y) (B.ll) 

J y°=const. 

$(x) = -i I d 3 y $(y)i°S(y-x;m) (B.12) 

J y°=const. 

with 

S(x; m) = {i$ x + m)A(x; m 2 ). (B.13) 
Furthermore we need the perturbative propagator 

(0| T4(x)$(y) |0) = Jd R p [6(x° -y )e- ip{x ~ y) (^+m) 

+6{y° - x )e +ip{ - x - y) (-tf+m)]. (B.14) 



C Calculation of M. 

The calculation of M. is demonstrated best in calculating the part 
/ = «0|T [d 4 x 3 [9(-x° 3 )A l ,(x 3 )e i ^ + ^ - d(x° 3 )A^x~ 3 )e i( ^ + ^} 

^^(-X?)^!)^!)?^!)-^^?)^!)^!)^!)] ... |0)) 



(C.l) 
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of the matrix element in (4.14). This expression is in principle governed by the four 



integrals I\, J 2 , I 3 and J 4 , which will be calculated successively in the following. We 
will start with I\ 



h = «0| T Jd^ei-x^A^x^+^e^A^) |0» 
= 0(-x°) [d%6(-x 3 )e l ^ + ^((0\ TA^x 3 )A p ( Xl ) |0» 



-ig„fl{-£) jd%0(-4)e^ J ^ Yk 2_ l X 2 + ie ^ k{X3 - Xl) 
-ig^i-xl) J d\ 5®(p[ +p 2 -k)J dx 3 e(-x 3 )e i( *° +p * 



dk° 1 

(2?r) k 2 -X 2 + ie' 



ik°(x°-x°)^ (C.2) 



The integration of k° is best done using the residue theorem. Unfortunately, because 
of the occurring ^-functions in x 3 and x\ the integral is neither in the upper nor in 
the lower complex plane convergent. This can be circumvented via 

d 4 x 3 6(-x° 3 )f(x 3 ) = [ d 4 x 3 f(x 3 ) - [ d 4 x 3 6(x° 3 )f(x 3 ) (C.3) 



so that we get two contributions, one in which the Xg-integration stretches over the 
whole real axis, we get a 5-function in k° and its integration becomes trivial and a 
second contribution where we can apply the residue theorem without any problems 

h = -ig w 9{-xl)[ l_ e <(w-n»>*i 

s — X 2 + le 



- [ d R k(27T) 3 5^(p 1 +p 2 - k)e ikxi - n =— 

J £>V + Vo - 



p\+p%-k° + ie 1 ' 

(C.4) 



With the help of 



J d 4 B5 {1) (B° - x° 3 )A(x 3 - B)choe~ ikB 



= e~ lkx ^ (e 2ik °' x °(k° + A; ') + (-k° + k Q ')e~ 2ik °'A with k ' = \Jk 2 + X 2 

(C.5) 

the integration of I 2 is straightforward. We can apply the residue theorem directly 
because the ^-functions ensure the convergence of the integral over k° 



«0| T J d%9{xl)A^x 3 )e i ^+^M-x° l )A p {x l ) |0» 
Jd 4 x 3 6(xly^ + ^m TA^xJApfa) |0» 
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= 0(-a;°) j dS^xlY^+P^ J d 4 B5^(B° - x° 3 ) 
A(x 3 - B)&«0| TA M (S)A p (xi) |0» 
= ig w e(-4) I d R k(27r) 3 5^(p 1+ p 2 - k)e^ Q 1 , . (C.6) 

To calculate / we need only the difference between both integrals 

[ 1— + * / ^(27r) 4 5( 4 )( Pl +p 2 - fc)] 

= -ig^i-xl) _ 1 . e 'fr+*»K (C.7) 

S /\ Id 

I 3 can be treated in the same way as I2 

h = «0| Tjd^ei-x^A^e^+^eix^Apix!) |0» 

= / d R k{2*r^\f l+ f 2 - ky kxi vUv * Iw-ie - (c - 8) 

The most complicated integral is h 

h = «0| Tjd%9(xl)A^£ 3 )e^ + ^M4)M^) |0» 
= jd^eixly^+^ity r TA ll {x^)Ap{x 1 ) |0» 

= 0(x°) Jd%9(x° 3 )e l( P^ x ~ 3 J d 4 B5^(B° - x° 3 )A(x~ 3 - B) 

d^CS^^C - x°)A(fi - C)((0| T&^(S)^oA p (C) |0)). (C.9) 



From this integral we get two contributions because of an additional contact term 
which arises due to the two time derivatives in the 2-point function 

«0| Td B oA,(B)d c0 A p (C) |0» = 

= d B od C o{{0\ TA^B)A p (C) |0» + ig w 5^\B - C). (CIO) 

The contact term gives rise to the contribution 

hi = -tg»AAV {pi+mYxi 

J d R Jb^L(27r)W(pl + p 2 - £)(e 4 ^? + e- 4ik °*° - 2) 

(C.11) 
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and the rest, again with the use of ( |C.3| ) and ( |C.5| ), results into 



_ e i(Pi+p 2 )x 1 ( /" J_(27r) 3 5( 3 )(pl + p 2 - jfe)(2 - e 4ifc0a; ? - e ~^*\) 

-)], (C.12) 



s — A 2 — ie' 



so that we get 
h = 



+e 



s — A 2 — ie 

Again we only need the difference between both integrals 

h-h = i gtip 6{ x i) ei(pl+P2)Sl 



(C.13) 



[- a_ +i / ^(27r) 4 ^( Pl +p 2 + 
f 



s — X 2 — ie 



(C14) 



so that we finally get 
f 



—2 



s — A 2 — ie 
T te 1 [^(-x?)g(x 1 )7 M g(x 1 )e J 



51 - ^(x?)g(x 1 ) 7 ^(£i)e l 



8 *]...|0». 
(C.15) 
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